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Abstract 

We completely solve the equivalence problem for Euler-Bernoulli equation using Lie 
symmetry analysis. We show that the quotient of the symmetry Lie algebra of the 
Bernoulli equation by the infinite-dimensional Lie algebra spanned by solution sym- 
metries is a representation of one of the following Lie algebras: 2A\, A\ © A2, SAi, 
or ^1.3,3 © A\. Each quotient symmetry Lie algebra determines an equivalence class 
of Euler-Bernoulli equations. Save for the generic case corresponding to arbitrary 
lineal mass density and flexural rigidity, we characterize the elements of each class 
by giving a determined set of differential equations satisfied by physical parameters 
(lineal mass density and flexural rigidity). For each class, we provide a simple repre- 
sentative and we explicitly construct transformations that maps a class member to 
its representative. The maximally symmetric class described by the four-dimensional 
quotient symmetry Lie algebra ^3^ © A\ corresponds to Euler-Bernoulli equations 
homeomorphic to the uniform one (constant lineal mass density and flexural rigidity) 
. We rigorously derive some non-trivial and non-uniform Euler-Bernoulli equations 
reducible to the uniform unit beam. Our models extend and emphasize the symme- 
try flavor of Gottlieb's iso-spectral beams (Proceedings of the Royal Society London 
A 413 (1987) 235-250) 
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1 Introduction 



Da Vinci and Galileo foresaw the need for a theory of vibrating thin beams. 
However they suggested theories that were either incomplete or erroneous. Da 
Vinci's theory was more descriptive and based on detailed sketches rather than 
physical laws and equations : he lacked tools such as Hooke's law, Newton's 
laws, and calculus which postdate him. In Galileo's approach, the nemesis was 
an incorrect calculation of the load carrying capacity of transversely loaded 
beams. We owe the first consistent thin beams theory to the Bernoullis. Jacob 
Bernoulli developed an elasticity theory in which the curvature of an elastic 
beam is proportional to its bending moment. Relying on his uncle elasticity 
theory, Daniel Bernoulli derived a partial differential equation governing the 
motion of a thin vibrating beam. Leonard Euler extended and applied the 
Benoullis theory to loaded beams. 

In Euler-Bernoulli beam theory, the transversal motion of an unloaded thin 
elastic beam is governed by the partial differential equation 



where f(x) > is the flexural rigidity, m(x) > is the lineal mass density, 
and u(t,x) is the transversal displacement at time t and position x from one 
end of the beam taken as origin. Equation (1) must be solved subject to initial 
and boundary conditions such as clamped ends, hinged ends, and free ends 
boundary conditions. 

In this paper, our focus is on the equivalence problem for Eq. (1): we seek 
necessary and sufficient conditions under which two equations of the form (1) 
can be mapped to each other using an invertible change of the dependent 
and independent variables. A particular case of this equivalence problem was 
tackled by Gottlieb [1] who was interested in equations of the form (1) that 
are equivalent to the uniform ( constant / and m) beam equation. In the same 
vain, Bluman and Kumei [2] (Chapter 6, Section 6.5) studied the problem of 
reducing a linear partial differential equation to a constant coefficient one. 

The layout of this paper is the following. There are four sections including 
this introduction. Section 2 deals with the complete Lie symmetry classifi- 
cation of Eq. (1). Section 3 is dedicated to the construction of equivalence 
transformations. We recapitulate our findings in Section 4. 




t > 0, < x < L, 



(1) 
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2 Symmetry analysis of Euler-Bernoulli equation 



Our goal in this section is to study the symmetry breaking of Eq. (1). We 
assume that the reader is familiar with the rudiments of Lie's symmetry theory 
[2,3,4]. 

A vector field 

d d d 

X = r(t, x, u) — + f (i, x, u) — + 7}(t, x, u) — (2) 



is a Lie symmetry of Eq. (1) if 



dx 2 



dx 2 



dt 2 



0, 



(3) 



Eq.(l) 



where X' 4 ' is the fourth prolongation of X which is calculated using the for- 
mulas 



_d_ 

du 



, J=(ji>h), \J\=ji+32, 



i<\j\<k 
uj = d lJl u/dt jl dx h , 

Vj = d j(v - TU t - Z,u x ) + r u J:t + ^ u J)X , u, lr = duj/ dr, 

Dj = D n D j2 , D n = (D t y\ D n = (D x )i, 

^ d d d d 

D t = ^- + u t — + u tt - h u tx - 1 

at ou ou t ou x 

d d d d 

D x — 7^ 1" u xT{ 1" u txT{ 1" u xx~ 1" 



dx 



du 



du t 



(4) 

(5) 
(6) 
(7) 

(8) 
(9) 



Since the symmetry coefficients r, £ and r do not depend explicitly on the 
derivatives of u, the left-hand side of Eq. (2) is a polynomial in the derivatives 
of u. Thus we may set its coefficients to zero to obtain an over-determined 
system of linear partial differential equations. In order to avoid the appearance 
of the integral /(to/ fY^dx in our calculations, we express the lineal mass 
density as follows. 

m(x) = (g'(x)) 4 f(x). (10) 



After some calculations, the determining equations for the symmetries simplify 
to 
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T = 4c 1 t + c 2 (11) 

? = ^ + T' (12) 
, = _(f^ + ^ + w + s^ },, +o( , (13) 

V /# fg 9 9 J 

(fa xx ) xx + g' A fa tt = (14) 

c 1 H 11 + c 3 H 12 = 0, (15) 

Cl H 21 + c 3 # 2 2 = 0, (16) 

( Cl H n + c 3 H 12 ) xx = 0, (17) 

where C\ to c 4 are integration constants, and the differential functions are 
relegated to Appendix A due to their size. It can be readily seen that Eq. (17) 
is a mere differential consequence of Eq. (15). Thus the equations we have to 
solve are Eqs. (15) and (16). 

It is well-known (see for example [2]) that the symmetry Lie algebra of a 
scalar linear partial differential equation is of the form U © where U is 
a finite-dimensional Lie algebra and L°° is an infinite-dimensional ideal of the 
symmetry Lie algebra spanned by the so-called solution symmetries. In our 

case 

L°° = < a(t,x)d u >, 

where a(t,x) solves Eq. (14) i.e. Euler-Bernoulli equation. Thus we are left 
with characterizing U . The determining equations contain four constants viz. 
Ci to c 4 with two constants that are unconstrained. Therefore the dimension 
of 17, r, lies between two and four. Below we elucidate all the possibilities. 



Case I: C\ and c 3 are arbitrary constants 



In this case, Eqs. (15)-(16) split into the following system. 

#ii = 0, (18) 

#12 = 0, (19) 

#21 = 0, (20) 

#22 = 0. (21) 

The system (18)-(21) is an over-determined system of nonlinear ordinary dif- 
ferential equations for / and g. Thus, a priori, we are not guaranteed of a 
solution. Replace Eq. (18) by the combination Eq.(18)-gxEq.(19). Solve the 
resulting equation for to obtain 

9 ~ 10 P 5 / + 2 g> [ZZ) 
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Use Eq. (22) to eliminate the derivative and from Eq. (19). It results 
the equation 

fW _/'/ (3) 1 11 r 2 i2/ /2 r 9 r r2S , 

7 / 10 / 5 / 2 10 / 3 1 J 



Employing Eqs (22)-(23) to get rid of the derivatives g^ to g^\ and 
from Eq. (20)-(21), we discover that Eqs.(20)-(21) are identically satisfied. To 
sum up, we have established that the over-determined system of Eqs.(18)-(21) 
is equivalent to the determined system formed by Eqs. (22)-(23). Thus pro- 
vided Eqs. (22)-(23) are satisfied, the finite part of the symmetry Lie algebra, 
L 4 , is spanned by the operators 



x *= u l> (24) 

d 2g d (f'g M Sgg"\ d 

X3=4t d-t + Y^-{w + ^) u ^l (25) 
X4 = ?^-l77 + ^J^ (26) 

Simple computations show that the nonzero commutators of the symmetry 
generators are 

[X 1 ,X 3 ]=4X 1 , [X 3 ,X 4 ] = -2X 4 . 
By making the change of basis 

e\ = Xi, €2 = —j^ X±, e 3 = - X 3 , e 4 = X 2 , 

it can be seen that L is equivalent to A 3j3 © A\ in Patera and Winternitz [5] 
classification scheme. 



Case II: C\ is arbitrary and c 3 = 



The determining equations (15)-(16) become 



#ii = 0, (27) 
#21 = 0. (28) 

We aim at rewriting the constraints (27)-(28) in terms of lowest possible deriva- 
tives of / and g. In order to achieve this goal, we first solve Eq. (27) for g^ 
to obtain 



5 



g W - 6 9" 9 {3) 2 9' 9 {3) 2 g' | fi <?" 3 | g g" 2 | 4 g" /* 



<7 ; g 5 / g' 2 g 5 / 

3 fv + g; f r _ 4 _ 3 + 3 . 

5 g f 5 gf 5 / 3 5 / 2 ^ 



(29) 



Employ Eq. (29) to get rid of the derivatives g^ to g^ from Eq. (28) and 
solve the resulting equation for /( 5 ) to obtain 



(5) -18 f 18 fV 54f 3 r _ 48/V/" _ 7/'/" 2 22<?'/" 2 



18 f A g" 48 f' 2 f'g" 22 f" 2 g" 22 f' 2 4f'g'f {3) 



5/V 5/V 5/#' 5/ 2 

16 /"/(3) | 2/7 (4) | (3Q) 



5/ fg' f 9 9' 

Thus the constraints Eqs. (27)- (28) are equivalent to the simplified constraints 
(29)-(30). Provided / and g fulfilled Eqs.(29)-(30), the finite part of the sym- 
metry Lie algebra, L 3,1 , is spanned by X±, X 2 , and X%. This Lie algebra 
corresponds to A 1 © A 2 in Patera and Winternitz [5] classification of lower- 
dimensional Lie algebras. 



Case III: c± — and c 3 is arbitrary 

Here, Eqs. (27)- (28) are equivalent to the system of equations 

H 2 i = 0, (31) 
H22 = 0. (32) 

By following the same modus operandi as in Case II, we arrive at the following 
simplified constraints on / and g. 



-3/'V , fg'f" 3f 2 g" , 4/V 



5/3 p 5/ 2 5/ 



9' 2 5/ 5 ' 
(5) = -18 f 5 54 f 3 f" _ 7f'f" 2 _ 18 fV 
1 • • / 1 5/ 3 P 5/V 

48f' 2 f"g" 22 f" 2 g" 22/' 2 / {3) 
+ 5/V 5/7 5/ 2 



+ , (33) 
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5/ fg' f g' { } 

The finite portion of the symmetry Lie algebra, L 3,2 , is spanned by Xi,X 2 
and X4. The Lie algebra L 3 ' 2 is nothing but the three-dimensional Abelian Lie 
algebra denoted by 3Ai in Patera and Winternitz [5] classification of lower- 
dimensional Lie algebras. 



Case IV: C\ = = c 3 



This is the generic case: there are no constraints on / and g. The finite part 
of the Lie algebra, L 2 , is generated by X 1 and X 2 . It is the two-dimensional 
Abelian Lie algebra 2A X . 

To sum up, we have established the following result. 

Theorem 1 Denote the symmetry Lie algebra of Euler- Bernoulli equation 
by S , and L°° the infinite- dimensional Lie algebra generated by the solution 
symmetries. Then, S/L°° = A^®Ai, Ai®A 2 , 3A 1; 2A 1 , depending on wether 
f and g respectively satisfy Eqs.(22)-(23), Eqs. (29)-(30), Eqs. (33)-(34), or 
are arbitrary. 



3 Equivalence classes and mapping to canonical elements 



From a symmetry standpoint there are essentially four classes of Euler-Bernoulli 
equations. Equations of the same class share the same symmetry structure, 
and are homeomorphic. Our aim in this section is to select a simple represen- 
tative for the non-generic classes (i.e. all the cases save the case where / and 
g are arbitrary), and constructively show how class members are mapped to 
their representative. 

In the sequel, we shall use capitalized variables to describe representative of 
each class. The rational behind this choice shall be apparent as this section 
unfold. We denote by [£] the set of all Euler-Bernoulli equations having C as 
the finite part of their symmetry Lie algebra. Notations not introduce here 
are those of the previous sections. 
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3. 1 The class [A 3>3 © A x ] 



3.1.1 Construction of similarity transformations 

It can be readily verified that f(x) = 1 and g(x) — x satisfy Eqs. (22)- (23). 
Thus, we select as representative element of this class the equation 

Uxxxx + U TT = 0. (35) 



The finite part of the symmetry Lie algebra of Eq. (35) is generated by the 
vectors 



The invertible transformation 

T = T(t,x,u), X = X(t,x,u), U = U(t,x,u) (37) 



maps an element of L4 3i3 ©Ai] to Eq. (35) if and only if the same transformation 
maps < Xi, X 2 , X 3 , X 4 > to < Y 1: Y 2 , F 3jl , Y 4 >. We look for a transformation 
(37) such that 

X^Yu X 2 ^Y 2 , F 3 ,i X 3 + HiXi + H2X4, Y A ^X A . (38) 



Recall that the transformation defined by Eq. (37) maps a vector field V 
depending on t, x and u to the vector field T(T) 8t + T(X) dx + T(U) du- 
Thus, in order to realize Eqs. (38a)- (38b), and Eq. (38c), we have to impose 



X 1 (T) = 1, X 1 (X) = 0, X 1 (U) = 0, (39) 

X 2 (T) = 0, X 2 (X) = 0, X 2 (U) = U (40) 

X 4 (T) = 0, X 4 (X) = 2, X 4 (U) = (41) 

Solving Eqs. (39)-(41), we obtain 



T = t + k u X = g + k 2 , U = huJfg'\ (42) 



where ki,k 2 , and k 3 7^ are constants. It can be readily verify that Eq. (38b) is 
satisfied for \i\ = k\ and ji 2 = k 2 . We summarize the findings of this subsection 
as follows. 
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Theorem 2 Equations o/L4 3j3 ©Ai] are homeomorphic to Eq. (35). A trans- 
formation that maps an arbitrary element of [A 3 3 © A ± ] to Eq. (35) is given 
by Eq. (42). 



3.1.2 Examples of non-uniform beams homeomorphic to uniform beams: gen- 
eralization and symmetry justification of Gottlieb 's iso-spectral models 

Here we look for closed-form solutions of the uncoupled system of Eqs. (22)- 
(23). 

We may rewrite Eq. (22) as 

to*} = lo 7* -5 7' (43) 



where {y,x} = y" jy' — (3/2) {y"/y') 2 is the so-called Schwartzian 'derivative' 
(it is not really a derivative but a differential invariant!) of y with respect to 
x. From the well-known result on inversion of the Schwartzian derivative [6], 
we infer that 

(44) 



where y\ and y 2 are two linearly independent solutions of the second-order 
linear ordinary differential equation 

= 0. (45) 




We are now left with solving Eq. (23). Its symmetry Lie algebra is spanned 
by the operators 

d d d 



Since the Lie algebra < ri,r 2 ,r 3 > is solvable (its second derived algebra is 
trivial), we may use successive reduction to depress the order of Eq. (23) by 
three to obtain an Abelian equation of the second kind which we could not 
solve in closed-form. For details about the successive reduction of Eq. (23), we 
refer the reader to Appendix B. Due to the lack of closed- formed solution of 
the reduced equation, we look for an invariant solution of Eq. (23). The most 
general invariant solution under a linear combination of Ti to T 3 are 

/ = K(Ax + B) m , f = Ce Dx (47) 
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where K, A to D, and m are constants. Note that the second ansatz Eq. 
(47b) satisfies Eq. (23) if and only if D — 0, and this possibility is already 
included in the first ansatz. The substitution of Eq. (47a) into Eq. (23) yields 
the constraint 



m(4m 3 - 32m 2 + 79m - 60) = 0. 



(48) 



Solving Eq. (48), we obtain 

f 3 5 1 
m e < 0, -, -, 4 > . 
I 2' 2' J 



(49) 



Using the ansatz Eq. (47a) into Eq. (46) yields 
„ A 2 m (4 - m) 



y" + 



20 (Arr + 5) 2 



y = 0. 



(50) 



The general solution of Eq. (50) is 



V = \ 



k 1 + k 2 (Ax + B) ifmG{0,4}, 

fci (Ax + Bfl A + A; 2 (Ax + Bf/ 4 if m e {3/2, 5/2}, 



(51) 



where fci and k 2 are arbitrary constants. We infer from Eq. (44) that 

if m G {0,4}, 



9 = S 



L+M(Ax+B) 
P+Q ( Ax+B) 



(52) 



L{Ax+B)V±+M{Ax+Bf/± r„ /2 ./ol 

( P(Ax+B) 1 /4+Q(Ax+S)3/4 11 171 fc l' 3 /^' /^/) 



where L, M, P and Q are constants satisfying LQ — MP ^ 0. 

The models defined by Eqs. (47a), (49) and (52) generalize Gottlieb's [1] iso- 
spectral models. 



3.2 The class [A x ® A 2 ] 



A particular solution of the system (29)-(30) that does not satisfy Eqs. (22)- 
(23) is f(x) = x, g(x) = x. Based on this particular solution, we choose as 
representative of [Ai © A 2 ] the equation 

(XU xx ) xx +XU TT = 0. (53) 
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The finite portion of the Lie symmetry algebra of Eq. (53) is generated by Yi, 
Y 2 , and Y 32 = 4T<9t + 2Xdx — Udu- An invertible transformation (37) maps 
a element of [A 1 © A 2 ] to Eq. (53) if and only if it maps < X x ,X 2l X 3 > to 
< Y 1: Y 2 ,Y 3>2 >■ We search for such a transformation by mapping the basis 
elements as follows. 

Xi-Yi, X 2 ^Y 2 , X 3 ^F 3i2 . (54) 

Simple calculations show that Eq. (54) is realized if and only if 

T = t + hg 2 , X = 2l 2 g, U^l 3 u^^, (55) 

where l±, l 2 ^ 0, and l 3 ^ are constants. 
Thus we have proved the following theorem. 

Theorem 3 Equations of [A 2 © A±] are homeomorphic to Eq. (53). A trans- 
formation that maps an arbitrary element of [A 2 © A{\ to Eq. (53) is given by 
Eq. (55). 

3.3 The class [3^] 

A simple solution of Eqs. (33)-(34) which does not satisfy Eqs. (22)-(23) is 
f\x) = 1, g(x) = \nx. To this solution corresponds the representative 

Uxxxx + X' 4 U TT = 0. (56) 

The quotient symmetry Lie algebra S/L°° of Eq. (56) is spanned by Y 1 , Y 2 , 
and Y 32 = 2 X dx + 3 U djj. An element of [3Ai] will be homeomorphic to Eq. 
(56) provided the Lie algebras < X 1 ,X 2 ,X i > and < Y ± ,Y 2 , Y 3j2 > are similar. 
We look for a similarity that transforms the basis vectors as follows. 

Xx^Yl, X 2 ^Y U X 4 ^Y 3 , 2 . (57) 

Elementary reckoning shows that the mapping (57) is realized by the trans- 
formation 

T = t + m 1 , X = m 2 e 9 , U = m 3 uy/Jg^9 , (58) 
where m\, m 2 ^ 0, and m 3 ^ are constants. 
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In summary, we have proved the following statement. 

Theorem 4 Equations of [3Ai] are homeomorphic to Eq. (56). A transfor- 
mation that maps an arbitrary element of [SAx] to Eq. (56) is given by Eq. 
(58). 



4 Conclusion 

We have studied in details symmetry breaking of Euler-Bernoulli equation. We 
have shown that the Lie symmetry algebra of the Euler-Bernoulli equation is 
one of the following: 2A 1 © 3Ai © L°°, A 1 © A 2 © L°°, or A 3t3 © A 1 © 
where L°° is the infinite-dimensional Lie algebra spanned by solution 
symmetries. Equations admitting a given symmetry class are characterized 
completely in terms of a determined set of non-linear ordinary differential 
equations that physical parameters (flexural rigidity and lineal mass density) 
must fulfill. Equations of the same class can be mapped to each other via 
invertible transformations. For each class we provided a simple representative 
and we explicitly constructed similarity mappings. For the particular class of 
equations equivalent to the uniform Euler-Bernoulli equation, we rigorously 
constructed explicit non-trivial examples that extend and generalize Gottlieb's 
[1] iso-spectral models. 
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Differential function appearing Eqs. (15)-(17) 



_6/' 2 6gf' 3 a , lOgf'f" 6gf' 2 g" 8gf"g" 30 / g" 2 

Fll -TT"7V " 8/ + T^ + ^^ + ^~ 

60/ffff" 3 4g/( 3 ) 20/g( 3 ) 60 f g g" g® 

-6/' 3 , 10 /'/" efV^S/V 60/fl" 3 4/W 60/</' g 
a 12 - ro , H ;— ; r- pr - H 75 ta ; r 



pg> fg> fg* g>* f g> g' 3 

io/g (4) 

_12f' 4 12 gf' 5 28 f 2 f" Ugf' 3 f" 10 f" 2 21 g f f" 2 
21 / 3 /V f 2 P9' f /V 

12 5 /' 4 5" , 6/' 3 g" 28 5 f 2 /" 5 " 1 1 /' /" g" 10 g f" 2 g" 



/V 2 /V /V 2 /<?' /ff' 2 

12 5 /' 3 5 " 2 6/V 2 22 g f f" g" 2 3 /" g" 2 12 gf 2 g" 



/V 3 iY 2 /<?' 3 </ 2 iY 4 

60 /V' 3 , 6gf" g" 3 , 120 5 /' </" 4 , 180 / g" 4 360 /<?<?" 5 



g' 3 + <?' 4 + g' 5 + </' 4 <?' 6 

10f/gj 12 gf 2 fW 9gf"f® 10ff/V/g 6 g" /( 3 ) 

+ 7 Tv + ~~ /v TP + 

12ff ff //2 /( 3 ) 6g/ 3 g( 3 ) 4/'V 3 > 11 gf'f'gW 2 f" g^ 
g' 3 + /V 2 /<?' /<7' 2 + 5' 

, 12 fffVffgj , 70 /' 5 " 5 ( 3 ) 6 5 /" 5 " g^ 180 g /' g" 2 g( 3 ) 



iY 3 g' 2 g' 3 g' 4 

300 /g"V 3 ) , 720 fgg»*g<® , 6g/( 3 )g( 3 ) , 30 gfg( 3) 



+ . + .o + 



</ 3 ff' 5 </ 2 <7' 3 



60 fg^r 270 fgg"gW 2 (4) Zgf'fW 4gg"/( 4 ) 

g' 2 g' 4 1 fg' g' 2 

2gf 2 g^ 15/Vj g/"gW 40 gfg"g( 4 ) 75 fg"gW 

fg' 2 g' g' 2 g' 3 g' 2 

180 /gg" 2 g (4) 6Qfgg<®gW g/ (5) 5gf'g {5) 12 f g {5) 



g' 4 + <?' 3 g' g' 2 

+ 30 /gg"g (5) 3 / g g( 6 ) 



g' 3 g' 2 

-12 f 5 34 f' 3 f" 21 f'f" 2 12 fg" 28 f' 2 f" g" 

22 ""^ + _ F? TV 7V~ + Pg' 2 

10 p 2 g" 12f 3 g" 2 22 f f g" 2 12 f' 2 g" 3 6 /" g" 3 
+ — - - + 



/<?' 2 Pg' 3 fg' 3 fg' 4 g' 4 
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120 /V 4 360 fg" 5 12/ /2 /( 3 ) 9/"/( 3 ) 10/V'/ (3) 

~l ^5 Tfi Jo J "I - 



„/6 



/V 



/<7 



/2 



12 5 //2 / (3) 6/ /3 g( 3 ) 11 f'f'gW 12 f' 2 g" g^ 6/V'<? {3) 



y3 



<r Pg rz fg' z fg' 3 

180 /V'V 3) 720 /g"V 3 ) 6 /( 3 ) g( 3 ) 30 f g^ 2 

~l" .9 ~l~ 



,/3 



^2 



/5 



V2 



V3 



270 fg"gW SffW Ag" 2/ /2 gW fVg 40 /W 4 ) 



.9 



.9 



,2 



g 



,2 



g 



,3 



180 fg" 2 g^ mfg^g^ /< 5 > 5 f g^ 30 f g" 3 f 

+ /3 .o + 



.9' 



/4 



g 



g' 



g 



,/2 



3 



/3 



.9' 



/2 



B Successive reduction of the order of Eq. (23) 

The Lie brackets of the symmetry operators are [ri,r 2 ] = r\, [ri,r 3 ] = 
0, [r 2 ,r 3 ] = 0. Thus < Ti > and < T 3 > are ideals of the symmetry Lie 
algebra. We may start reduction using Ti or T 3 . It is crucial start reduction 
using an ideal of the symmetry Lie algebra in order to preserve the remaining 
symmetries. 

A basis of first-order differential invariants of Ti is formed by / and /'. We 
define new dependent and independent variables by 

y = f, t = f. (B.i) 

In the new variables, Eq. (23) reads 

y ~ t y 5 t 2 + 20 y y 2 + 10t 3 ' [ ] 

where the over-dot stands for differentiation with respect to t. In the new 
variables, the symmetries T 2 and T 3 are 

T 2 = -yd y , T 3 = td t -ydy. 

It can be verify that Eq. (B.2) inherits the symmetries T 2 and T 3 as expected. 

A basis of first-order differential invariant of T 2 is formed by t and 

z = V- ■ (B.3) 

y 
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In the new variables t and z, Eq. (B.2) becomes 
•■ z n . 12 41 z 2 ir 3 9 1 

z = 7zz z^ 16^ 3 H -• (B.4) 

t 5 20 t 10 t 3 K ' 

In the variables t and z, the symmetry operator T 3 becomes 

T 3 = td t - z d z . 
A basis of first-order differential invariant of T 3 is 

u = tz, v = tz + t 2 z'. (B.5) 



In the new variables u and v, Eq. (B.5) is 

dv „ „ 320m 3 -181m 2 + 108m-18 

= 5 _ 7 M . (B.6) 

du 20 v v ; 



Eq. (B.6) is an Abelian equation of the second kind. 
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